# Full text of "Microlensing in phase space II: Correlations analysis"

## See other formats

Mon. Not. R. Astron. Soc. OOO.H1IT7IC2006') Printed 5 February 2008 (MN I*TeX style file v2.2) Microlensing in phase space II: Correlations analysis A.V. Tuntsov* & G.F. Lewis f A29, School of Physics, University of Sydney, NSW 2006, Australia \0' o o <: > O o o Oh 6 % Accepted 2006 April 6. Received 2006 April 4; in original form 2006 January 16 ABSTRACT Applications of the phase space approach to the calculation of the microlensing auto- correlation function are presented. The continuous propagation equation for a random star field with a Gaussian velocity distribution is solved in the leading non-trivial approximation using the perturbation technique. It is shown that microlensing modu- lations can be important in the interpretation of optical and shorter-wavelength light curves of pulsars, power spectra of active galactic nuclei and coherence estimates for quasi-periodic oscillations of dwarf novae and low-mass X-ray binaries. Extra scatter in the brightness of type la supernovae due to gravitational microlensing is shown to be of order up to 0.2™ depending on the extent of the light curves. Key words: gravitational lensing - methods : statistical stars : oscillations - galaxies : active - dark matter methods : analytical 1 INTRODUCTION Statistical approaches to gravitational microlensing have a long his tory th at extends back to at lea st the mid 1980s (e.g.JPaczvnski 1986, Kavser. Refsdal fc Stabell IQSl IPeeuchi fc Watsonlll987fl when it was realized that in gen- eral, a large number of stars contribute to the deflection field, and as a result, it is nearly impossible to model in- dividual micro lensing light curves in detail unambiguously JSchneider fc Weiss 1987: Katz. Balbus fc Paczvhski 1986) . Since then, a number of methods have been developed to cal- culate the statistical properties of microlensing light curves. Apart from a few early p apers, the methods employed are numerical - ray-tracing ([K avser, Refsdal fc StabcU 198d eisj Wambseanss. Witt fc Schneider. .1992 : Schneider fc Weisd 1983) and image-tracking techniques JLewis et alJ Il993 i Witdll99a) - because they alone are capable of calculat- ing the magnification probability distribution function over a large range of magnification values. However, it was soon realized that analytic methods could help solve the complementary problem of correla- tions in the magnification values at nearby positions in a microlensing ligh t curv e. Based in part on the work of [Peguchi fc WatsonI lll987l) . ISeitz fc Schneided lll994) deveb oped a method to calculate the deflection angle probability distribution function and successfully applied it in a deter- mination of the flux autocorrelation function for a source observed through a single lens plane for a variety of the pa- rameters iSeitz. Wambseanss fc SchneideJll994l) . However, this method did not attract much further interest because of * E-mail: tyomichOphysics. usyd.edu.au t E-mail: gfl@pfiysics.usyd.edu.au substantial technical difficulties in its implementation. Nein- dorf (2003) took a further major step forward by developing a simple approximation to the correlated deflection proba- bility distribution that could be applied to many interesting astrophysical situations. In a companion paper (Tuntsov fc Lewis 2006, hereafter Paper I), we presented a similar approach in which a phase space (or intensity) description of the deflection process is used before projecting intensities onto observable quantities. The method makes clear some of the important approxima- tions used in studies of gravitational lensing and is easily generalized to include many lens planes along the line of sight, which is especially interesting in the case of a cosmo- logically distributed population of lenses. To avoid arbitrary 'slicing' of the space into individual lens planes, we have also proposed a continuous approximation to the propagation of intensity moments arguing that for average quantities there is no conflict between the thin lens approximation and a 3D distribution of lenses; we call this 'nearly' 3D lensing. In Paper I we have also solved the problem of the prop- agation of flrst-order moments. However, this is of little in- terest by itself as the results follow trivially from the conser- vation laws and the required symmetries of the initial condi- tions. In the present paper we apply these results to the non- trivial problem of the propagation of second-order moments of the flux; traditionally, the most interesting projections are the autocorrelation function and its Fourier transform - the power spectrum. In section 2, we review the results of Paper I in relation to this subject and present a general model for the initial conditions. In section 3, the deflection potential is calculated in the flrst non-trivial approximation for a random star field with a Gaussian distribution of velocities. The formal so- Tuntsov & Lewis lution of the propagation equation using the perturbation technique is obtained in section 4, which is followed by the calculation of the flux power spectrum for a point- like source in the leading order of the deflector density. For clarity, most of the calculations of this section are presented in the Ap- pendices. Section 5 discusses the application of these results. First we consider modifications to the power spectra of periodic and quasi-periodic sources. We show that known microlens- ing populations of the Galaxy can introduce ~ O.f^-deep modulations in the pulse profiles of a kpc-distant pulsar and significantly affect coherence estimates for quasi-periodic os- cillations of X-ray binaries and dwarf novae in the mHz to Hz frequency range. Then, corrections to the autocorrelation functions of aperiodic sources are considered; in particular, we find that in the 0.01 — 0.1 range of redshifts used for the calibration of type la supernovae brightness, the extra scat- ter due to gravitational lensing can amount to 0.04™ — 0.2™ depending on the light curve sampling time-scales. We con- clude with a short list of refinements to the method that can increase the quality of our predictions. 2 EQUATION FOR THE SECOND ORDER MOMENTS In paper I we have considered the transformation of the statistical moments of the intensity field as light propa- gates through a universe filled with gravitational lenses. The second-order moments are defined as the average value of the product of intensities at a pair of positions in the phase space 7V/(^ 9 (ri, r2, ai, 0:2, ti,i2), defined by the spatial po- sition of a light ray in the lensing plane r^, its direction a.i and the moment of observation ti, {i = 1,2): M2 (ri,r2,ai,Q;2,ti,i2) = E J (ri, ai,ii) J (r2, 02, ^2) • (1) The expectation is taken with respect to an ensemble of realization of the intensity field; further discussion of this ensemble is given below. It was shown in Paper I, that the conservation of surface brightness in gravitational lensing - that is, the conserva- tion of the intensity along the ray - results in the linearity of transformation on M2 as the light rays propagate from one lens plane to another and get defiected at them. It is therefore convenient to consider M2 as a vector in some lin- ear function space V (M^\. In particular, M2 can be repre- sented in any other functional basis; certain properties of the transformations often suggest particularly convenient bases. Under the thin lens approximation it is possible to in- troduce two operators, termed deflection "I? and focusing F, which describe how the intensity moments are modified when the light rays are broken on the lens planes ($) and when they propagate freely between the planes {F). Then, we derived a Fokker-Planck-type equation to describe the continuous limit of this treatment, and to answer the ques- tion of how the intensity moments behave as a function of depth D along the line of sight: dM2 dD d$ Id dF dD Here d<i>/d_D and d-F/d_D are effectively the generators of <!> and F, describing the changes in M2 due to these operators applied over an infinitesimally small distance dD. In the case of the second-order moments, the focusing operator in the above equation remains the same: dF _ _ >r^ df5 ^^ 2^ =1,2 Oi dr + 2^a, (3) Unlike the first-order case, there is a difference in the symmetry properties of the deflection operator. The thin lens and small angles approximations still apply, and there- fore the deflection probability depends only on the deflection angle /3i = a'^ — cti, and not on a;; this leaves the associated deflection operator diagonal (proportional to a (5-function) in Ti coordinates - Fourier conjugates to cti. In the focusing operator, we simply replace cti with —id-n- At the same time, deflection probabilities are no longer entirely independent of temporal or spatial coordinates as the deflection angles and potential time delays of the two rays are correlated. Dependence on t means that the asso- ciated defiection operator is no longer diagonal in the tem- poral frequency basis. However, recalling that the potential time delays turned to be unimportant in the observationally accessible range of frequencies, we will ignore them in the subsequent analysis and in doing so avoid a difficult integral term in the equation for M2 . Then, the defiection probability operator is diagonal in ri,Ti,ti basis. In the stationary and spatially uniform case, probability densities for the light rays to be deflected by (3i depend on the difference between the spatial and temporal coordinates associated with the deflections: <l(ri,Ti,ti,r2,r2,i2 -> ri, ri,ii,r'2, T2, ^3) = (4) 'I>(ri,r2,ti,i2,Ti,r2)(5(Ti - ti)S(t2 - T2) x5{r', - ri)5(r'2 - r2)3{t[ - U)5(t'2 - ^2), so the equation for the second moment assumes the following form: dDM2 (D;ri,r2,ii,t2,Ti,T2) = / 1 ' 1 (5) 1 + TT^ii^t; 2^. .1^.1 ' 2c "^*=^ +P2 (ri,r2,ti,i2,Ti,T2) JM2 (_D;ri,r2,ti,t2,Ti,r2) where p2 is the derivative of $ with respect to D. To complete the statement of the problem, one also needs to specify the initial conditions M2|d=o = ^2- In the present study, we will restrict ourselves to the case con- sidered in Paper I: an isotropic source located at D = —Ds- Then, for the first moment one has: Mi'(k,a;,r)=/(k,a;)— -exp ujDs '"^ (r + kD,)^ 2luD. (6) M2. (2) where f{\i.,u)) is the Fourier transform of the source profile. However, there is an important difference between the first and second order moments. In Paper I the ensemble over which one should average to obtain the moments was not fully defined because it served as a conceptual tool in the general derivation. Two alternatives were provided for the definition - one either looks at independent patches of the sky and averages quantities of interest along different lines of sight or observes a certain source long enough to consider Microlensing in phase space II: Correlations analysis 3 different patches of the Ught curve as independent elements of the ensemble. In this paper, we are more concerned with practical ap- plications and for this reason we will stick to the latter defi- nition as it is closer to the situation encountered in observa- tions. In this case, an internally consistent definition of the averages can only be provided for sources that are statisti- cally stationary. In the first order, this leads to the trivial result - in the absence of a special point in time, all modes with (jj 7^ have zero amplitude and stay so because of the homogeneity of equations, while the u — mode is preserved due to conservation of energy; even for the alternative defi- nition using different sky patches, it is hard to think of any reason for a certain time position to stand out, as noticed in Paper I. This can be expressed in a formal way. Indeed, shifting the origin of the time axis by a constant T results in a phase factor e'""^ in front of the Fourier transform of a function F{t): iujT , Ft{uj) FH; (7) the same applies to all averages. Therefore, in a statistically stationary case, when EtF — EF for all T one has ErFiuj) = e'""EF(a;) =EF{iu), VT. Hence, EF{uj) = for all tj / 0. Similarly, for the second order moments one has ETF2iLOi,u;2) = e'^'^'+'^^^^EF2{uji,iU2), (8) (9) and the only non-trivial statistically stationary components are those with a;2 ~ — i^i = i^/2; we presume w ^ from now on. When F{t) itself is a real quantity - as is the case for intensity and flux - F{—lj) = F*{uj), and the function Pf{uj) = EF(-lj/2)F{lu/2) = E \F{Lu/2)f (10) is the power spectrum of F oscillations and represents a Fourier transform of the autocorrelation function of F. If one now assumes that the source brightens and dims synchronously over its entire spatial profile so that it can be separated into a temporal and spatial parts: /(k,<.^) = <7(k)/iH (11) the initial conditions for the second moment of intensity distribution of an isotropic stationary source located at _D = —Ds can be written as: M2(kl,2,fi,C^,Ti,; ■s{n) Ac^P{iu) X exp ljDs ((n + kiD, 5(ki)ff(k2) (12) -{T2+k2Dsf) where Q — {uji+uj2)/2, lo — lu2—uji. The average g{'ki)g{'k2) is used because correlations may be present in the source profile. From now on we will drop the factor 5(0.) which arises from the infinite extent of time axis and thus can be considered as a constant through regularization. 3 DEFLECTION POTENTIAL In calculating the deflection potential p2 we will again assume that the distribution of the deflectors is statisti- cally uniform and stationary and ignore any correlations between individual deflectors; we will also assume Poisso- nian distribution for the number of deflectors in a lensing plane. In t his case, the Markov-Holtsmark-Chandrasekhar iChandrasekharnl943i : jNeindort, 20031 method gives the fol- lowing result: P2 (Ti,T2,r,t) = iy\\C\\ (p(Ti,T2,r,t) - 1) , (13) where f is the spatial density of the deflectors, r = r2 — ri , t = t2 — ti, \\£\\ is the area of the lensing plane, and p is the Fourier transform of the defiection probability density with respect to the two deflection angles, calculated for a single deflector in the plane. The probability density itself can be easily obtained since, for a single deflector of mass m, the deflection an- gles uniquely determine its positions pi,2 with respect to the rays: pi — Vi— mfii/ 01. Therefore, for a uniform plane, ™4 p(/3i,/32,r,i) = /3!/3t \C\ (j) \ m /3? ||-^)+r,i)(14) Here the fraction is the Jacobian of the pi,2 -^ /3i,2 trans- form, while <j!>(p, t) is the probability density for the position change p over a time interval t - essentially, the velocity distribution of the deflectors. The simplest choice will be to assume a Gaussian distribution of these velocities with average v and dispersion g: r /_ ..-in2" (p,t) = Vl-Ko t j exp (p-vty 2aH^ (15) However, Fourier transforming 11411 even in this sim- plest case is a very hard task. We can, however, expand it in Taylor series in ti,2 and restrict ourselves to the most im- portant contribution in this expansion. This is a necessary exercise given that the problem for l|^ is in ten dimensions - numerical methods will not be easy to implement here unless we have some analytic understanding of the solutions. From the normalization of the probability density, the zeroth-order term in Taylor expansion of p(ti , T2 , r, t) is unity and it cancels out in 11311 . The coefficients in front of Ti and T2 are the average defiections of the two rays which equal zero for a uniform lensing plane. Therefore, the first term of the expansion is of the second order: p(Ti,T2,r,i) '--2 {{{t, ■ f3^f) + {{t2 ■ 02)')) (16) +2((Ti-/3i)(T2-/32))(r,i) + o{ti,T2] where the angled brackets correspond to averaging over the distribution 11411 . For symmetry reasons, the averages on the flrst line do not depend on r and t, and therefore they cannot depend on the direction of ri,2 either. Thus, (((n . /3i)^) + {{t2 ■ (32)'}) = i{/3^) (rf + r^) . (17) The other average in H16^ is responsible for the correlations. It can be rewritten in a simpler form 2{(ti ■/3i)(t2 •/32))(r,i) = {P'}t^Co ir,t)T2, (18) where Co is the correlation matrix C;Mr,t) = 2{/3^)-^{/3i/3^)(r,i) (19) The average in the angled brackets can be found using the joint distribution function of /3i,2: Tuntsov & Lewis {PiPi)ir,t) d2/3id2/32p(/3i,/32,r,i)/3l/3^ (20) (1) Later we will need the the Fourier transform of the correla- tion matrix as a function of r and t: 2 ,, exp [i (k ■ r + cjt)] Co(k,a;) (27r)3 '-Co{r,t) (21) For a Gaussian distribution of velocities HSU this integral can be done analytically. One first integrates with respect to r (or equivalently, r — ?nAl//3) over an infinite plane to obtain a Gaussian-type integral with respect to t, while the integrals with respect to /3i,2 decouple. As a result, one obtains a degenerate matrix: where a marks the direction of k: k = kea and 2 m^ C(k,w) = ^/---^exp (v-k- 2cr2fc2 (22) (23) x[Jo(pofc/2)-Jo(Apofc/2)]^ where we have assumed that the upper limit m/po on (3 is determined by the lens physical size po while the lower one, po/R comes from the extent of the lensing plane R — Apo- For relativistic compact objects, such as black holes and neutron stars, po/m ~ 0(1) while for white dwarfs po/m ~ O(IO^). In the case of continuous distribution of the deflecting surface density E(p, f), the deflection angle power spectrum can be related directly to the statistical properties of E. For a spatially and temporally uniform case one has (E is in mass per length squared units): (/?i/?^)(k,^) 4G^ fc'fe^ PE(k,(.^) (24) C* fc4 where Pe is the projected density power spectrum, which can be found using an appropriate form of the Limber's equation. In the formula above we have neglected retarda- tion effects, therefore it is valid when the motion of deflecting matter is slow compared to the speed of light. However, the smooth-matter description requires care- ful treatment of the correlations along the line of sight and will not be considered in this paper. 4 FORMAL SOLUTION AND THE FLUX POWER SPECTRUM Neglecting the higher-order terms in Hlfcil , we may write JSJ as 9dM2 l/:||(/3')TfCo(r,t)T2Af2 (25) + ( 19x1 ■ 9ri + 19x2 ■ 9r2 + TT^^l ^«1 + TT^l^ ^t^ ] ^2 2c 2c -j\\C\\{f3'){T?+Ti)M2 In the absence of correlations represented by the first line of the right-hand side a dense system of eigenvectors Ma (and corresponding eigenvalues \a) of the remaining operator is provided by various products and associated sums of the form M4ri,r2,Ti,T2,ti,i2) = MW(ri,Ti,ti)Mi^'(r2,T2,i2) (26) Aa — y^ai + ^a xM (27) with a composite index a = (01,02), where MaJ , XaJ are the solutions for the first-order equation: 1 X'a'/M^l> (r, T, t) = [id^ ■ dr + y/^dtj M^J' (r, r, t) (28) -Jl|/:|l(/3^)r^MW(r,r,t). As was found in Paper I, these represent harmonic modes in r and t multiplied by an angular functions Cmni''') depending on the spatial and temporal frequencies k, uj as parameters (see Appendix IXt. The functions Ma^ form a complete orthonormal ba- sis of the solutions, while their products of form 12611 do the same job for the problem considered in this paper. It is therefore convenient to rewrite 12511 in the basis of Ma- Namely, let M2(D;ri,2,Ti,2,ii,2) = Vcae^°-°M4ri,2,Ti,2,ii,2) (29) with a — (ki,2,<jJi,2, ?wi,2,'T.i,2) running through all possible values. When there are no correlations a set of Ca — const chosen such as to represent the initial conditions at D = gives a solution of 12511 . In the presence of correlations, which we will treat as a perturbation, we allow these 'constant coefficients' to vary with D to obtain the following system of linear equations in the perturbed case: Ca — / ^aa'^ (A ,-\a)D (30) where the dot denotes differentiation with respect to D. The matrix elements Caa' are defined as Caa' = -M\{(3'} {MlT^CoT2Ma') (31) The perturbation method then proceeds by introducing a dummy parameter e to characterize the order of perturba- tion eC and expanding Ca{D) in its power series: Ca{D) = c°(Z?) + €cl{D) + e\l{D) + . (32) such that (^{D) functions correspond to no correlations at all, c\ (D) are proportional to the first power of perturbation and so on. A sequence of successive corrections to Ca{D) can now be obtained by equating coefficients at common powers of e: Ca = 0, C„„/e' (x^,-\a)D n (33) (34) and so on. From the first of the above equations cJJ = const and if we choose these constants to represent the initial conditions M2 at D = 0: c° = (m],M2), the initial conditions in aU other orders will be uniform 4(0) = 0, i€jV. (35) Physically, the zeroth-order term describes just the partial loss of the phase coherence in different modes as a result of geometric time delays due to uncorrelated random deffec- tions, as explained in Paper I. Microlensing in phase space II: Correlations analysis 5 Further, we restrict ourselves to the first non-trivial term in the expansion 13211 . It is proportional to the strength of correlations in the deflection angle at different points and describes how these correlations translate across to the in- tensity field. By mixing the intensity modes of one- ray statis- tics CaiCaJ exp[(Aai -|- Aaa)-^] in a correlated way, extra cor- relations in the combined field are created, which can be described by solving 1341 using c° = const: e c„(D) = > c„, dD' Ca'e XaiD~D') X iD' (36) The power in the first-order term is also turned into extra signal by correlated deflections, but the higher-order contri- butions will be neglected in this study. Although the intensity field provides a complete de- scription in the sense discussed above, the intensity is rarely observed in the real life due to the limited angular resolution of our telescopes. The quantity that is normally measured is the fiux F. _F is a linear functional on the intensity field and, just like intensity, defined for every spatial and time coordinate r, i. Thus, in calculating it, one only needs to consider the projection in the angular sector - namely, consider t — 0. We can therefore prescind from the exact parametrization of r and t dependence - e.g., the temporal Fourier harmonics of the flux observed at points ri,2 -F2(ri,r2, 1.^1,(^2) = (37) will be those of the intensity projected at t = without any additional transformation. In the approximation we use, there are two contribu- tions to these quantities, associated with the zeroth and first-order terms in the perturbation theory expansion. They both depend linearly on the initial conditions and there- fore one can introduce two D-dependent functionals Fj ' £ V''{M^) that calculate the flux Fourier components as func- tions of M2 at respective accuracy. In the Appendix]^ these functionals are calculated with respect to Ma functions. Since they form an eigenbasis of the unperturbed problem, F° is diagonal in this representation. This is not the case for F^: correlated deflections can pro- duce the power in the flux Fourier spectrum at frequencies where the flux power spectrum due to initial conditions is zero. A stationary, isotropic and synchronous source is de- fined by just two functions: the average source spatial profile Fourier transform g(ki)g(k2) and fluctuations power spec- trum P{u)). Therefore the expectation value of the autocor- relation function Fourier transform (ki,2 = K =F k/2): F°'^ (K, k, u^;D)=Y. i^:,}(23)P(a;')<?(k'i)<;(k'2) (38) with d — (K,k,aj); here we use the summation instead of the integration sign to avoid unnecessarily long notation. The coefficients F^^, (D) are given in the Appendix ^ Eis an explicit expression for F^^, and as an integral of an explicit expression for F^^,. They describe the fraction of variability power that, in the presence of correlated defiec- tions on path length of D, is transferred from the original (a;', K', k')-mode to a mode at (cj,K,k). As noted above, F^ is diagonal in this basis. However, a question that is more interesting from the observational point of view concerns the value of the fiux temporal variability power JT observed at a certain point R (or autocorrelation between a pair of points ri_2 = R=F i'/2) in the plane of the observer. This can be obtained by sum- mation of all amplitudes multiplied by their Fourier mode value at the given point (s). ^"•^(^,0)= / dLj'd'kid'k2d^k'id k2F^'^'(D) (39) xg(k;)g(kyp(c^') ,, exp[-i(ki • ri) - i(k2 • r2)] 47r2 For the non-trivial term, coefficients F^^, depend not only on ij,ui' but spatial wave numbers as well, and therefore the summation can only be done for a certain source profile For the purpose of illustrating this method we choose the simplest case of a point-like source at the origin of the source plane and a single observer at the origin of the plane of the observer ri,2 = 0. To facilitate comparison with the results of Paper I, we assume that the source is located at D — —Ds ^ and that D > space is filled with deflectors of constant number density u. We perform the calculations in Appendix |n| here only the necessary results are present. In order to visualize the results we consider a 'line pro- file broadening' function ^{lj,uj' ,D) that represents broad- ening of a single line in the power spectrum at a certain u': P{cl>) = P{iu')5{iL!—uj'). 4> is defined as the ratio between 1391 and the value of the flux power spectrum T^{u',D)\ which would have been observed at D were the lensing pop- ulation absent. The latter can be obtained by noticing that the power spectrum is quadratic in flux which obeys the inverse- squares law when the light rays are not deflected. For the above choice of the source profile and the position of the observer (see the Appendix!^ .F°(^',g)^ ./^""l. ^0(v^^ {Ds + DY (40) and therefore i^°'^{Lo,uj',D) = {Ds + Df /d'Kd^K'd'kd'k'i?°^J(D)(41) The part of JT" due to uncorrelated deflections is of the second order in the deflectors' density v. As shown in Pa- per I, in all relevant astrophysical situations the density of deflectors is very low in terms of a dimensionless quantity rnD^JvLo/c that controls the strength of gravitational lens- ing by a three-dimensional distribution of compact lenses. The dominant non-trivial contribution to ^j) in this regime is supplied by the term due to correlated deflections J-^ ~ Olyy/v). Assuming that parameters u,m,a,v are con- stant along the line of sight one can introduce the frequency difference scale Aojo 2(j Po' (42) 6 Tuntsov & Lewis an amplitude parameter 3c Po (43) m 8n\nAv {D^+Dy and write down the following approximation for i/) (we drop superscript '1' because this is a dominant non-trivial contri- bution) : 12 exp |"-«V2cr2] 7/;(Ac^, D) ^ ^ i * ' Ds /27rlnAAa;o KAuJo'Ds+D'a where a scaled version of the line profile function -,g (44) *(i/,s,^,gj q I dxw {s,xq) exp y '2a:2 (45) Io(-^)[Jo{x)-Jo{Ax)]\ \(J X J An auxiliary function it)(s,^) is defined in the Appendix IHl and proportional to ^~^ for ^ 2> 1. Since A 2> 1 we can safely neglect the second Bessel function in the square brackets when evaluating (1451 except for some extremely low values of y when it does induce a small correction; Jo can also be neglected in subsequent formulae 15UII and 1521 . Figure displays the behaviour of the line profile ^ as a function of the frequency difference variable y for s = and w/o" = 1. It has a flat central part and then falls off gradually. The behaviour of '^ is rather hard to study ana- lytically as it has an essential singularity in the integrand; nevertheless, it can be understood qualitatively by noting that exp[— y^/2a;^] effectively cuts out x < y from the in- tegration region and using known form of w(^) and Bessel functions. There are two different regimes for <&(?/). At s = 0, for q < 0(1) the height of the flat part is proportional to q^ and it extends to y ~ q~^ followed by a '^{y) ~ ci~^y~^ tail afterwards: *(y) q V y<0(q~^) y>0(q-') q<0{i) (46) For q > 0(1) the height of the central part is oc g and the extent is again 0{q~^); however, the fall-off rate here is up to J/ ~ C(l) where it turns to the same *(y) <? ^y~ *(y) ■^ q y ^ tafl: y<oiq-') z y t >o(i) 0{q-^)<y<0{l) , g<0(l).(47) y . Oscillations of the Bessel functions enter the region not dom- inated by the exponential fall-off at different values of y which results in a few 'wiggles' in the resulting 'i!{y) curve seen in Figure In the opposite case of a distant source viewed through a localized overdensity of lenses, 1 — s ^ 1, the behaviour is generally the same but q parameter should be replaced with its scaled version q =q{l-- s ,3/2 m / Stt In A vD^ po (48) Figure 1. The line profile function ^{y, s,v/a,q) as a function of frequency argument y for s = 0, v/cr = 1 and a number of different values of the amplitude parameter q. The profile has a fiat central part of the height oc q (for q > 0{1)) or q^ (otherwise) extending to yc ^ 0{q up to y ~ C(l) and as q~^y~^ afterwards. Certain 'wiggles' can be seen in the shape of ^ which stem from the behaviour of the Bessel function in 1451 . the scaling arises from the fact that the main contribution to w(s,^) comes, for ^ ^ 1, from the region of i such that (l-t)^(t^ <r Different cuts through 'i/{y,q) can be obtained by con- sidering ^I" as a function of q, as shown in Figure |5| here, again, we assume that the source is inside the deflectors- filled volume: s = and choose v/a = 1. For a given fre- quency difference parameter y, <!/ first increases {y is a part of the central plateau) as either <& ex g^ for q < 0(1) or ^ oc g for g > 0(1), and then decreases as q~^ . The enve- lope curve of the ^(g) family thus represents the height of the central part as a function of the amplitude parameter g. Taking into account the scaling (1481 . s 7^ (to be more precise, 1 — s ^ 1) curves can be obtained by a shift of the ones given in Figure |5| to the right. For completeness, let us also state the dependence of the line profile on the velocity parameter v/a. At y <ti 1 ^ is nearly independent of v/a while at the shoulder and tail it scales as roughly exp[v'^ /2a^]. As a result, the shoul- der and tail of the product exp[—v^/2a^]'^{y,s,v/a,q) are insensitive to the velocity while the centre of the profile is significantly suppressed for v/a 3> 1 and gently connected to the shoulder by a practically featureless curve. The total fraction of the initial power added to the spec- trum by correlated deflections can be obtained by integra- tion of tp with respect to Alu: y{D)= / dAujiP{ALu,D) (49) 12exp[-«V4a2] ^ f v^ \ ,, f i^A ^°(,4^j^l Ds Ds + D' where 00 Y{s, q) = Jd^ ^w {s, [Jo i^/q) - Jo (AC/g)]' (50) Microlensing in phase space II: Correlations analysis 7 - - Y(q) ., Y,(q) Y,(q);q ,,, Y^(q)/qY(q) ^-^ 0.01 0.1 1 10 100 1000 10000 Figure 2. The line profile function ^{y, s,v/a,q) as a function of the amplitude parameter q for s = 0, v/a = 1 and a number of different values of y. For every frequency argument value y, function ^(q) reaches a maximum at qdy)'- Vciqc) = J/ and then falls off as y enters the wing of the line profile. The envelope curve of this family represents the value of * near the edge of the central part ydq)- Figure 3. Functions Y{s,q), Yi{s,q), Yi{s,q)/q and the dimen- sionlcss effective 'line width' Yi (s, q)/qY{s, q) as a function of the amplitude parameter q for the source embedded in the deflectors- filled volume: s = 0. Y and Yi both saturate for large values g 3> 1 and have oc q tails for g <§; 1. The dimensionless effec- tive 'line width', on the contrary, is a monotonically decreasing function of q. Another quantity of interest is the effective 'line width' defined as: (Auj) (D) = y'^ dAiu \Auj\ i/. {Aiu, D) where (51) ^ /T exp [^V4a^] y,(s,g) '^"V2 hiv^/^o^) qY{s,q) Y,{s,q) = / dUw{s,0 [Jo (C/g) - Jo (AC/g)] = (52) we need to make use of the absolute value of Alu in such a definition because as can be seen from the asymptotics of w{s, ^) and Bessel function the second moment of Aui does not exist. Using II45I I one finds that, for q < 0(1), the total frac- tion of the power created by correlated deflections y is pro- portional to q and is divided in roughly equal parts between the central plateau and y~^ tail. They also make approxi- mately equal contributions to the integral with respect to y in 1511 1. which is independent of q. As a result, the 'effective line width' (Ao;) oc q~^ . For q > 0(1), we find, using 1471 . that y saturates for g 3> 1; most of the power is split roughly equally between the central plateau that is high (oc q) and narrow (oc q~^), and the y~^ shoulder that extends out to Alo ~ 0(A(x;o), only a small fraction ~ 0{q~^) is in the y~'^ tail. The 'effective line width' in this case (Aw) oc In q/q. Figure |3 shows the results of numerical integration in (I50II52II . They agree well with the approximations above. An important property of the broadening function i/; in the low-density approximation is that it only depends on the frequency difference Au = u' — uo. The convolution the- orem then makes the transformation of the autocorrelation function (which, by virtue of the Wiener-Khinchine theorem is the Fourier transform of the power spectrum) very sim- ple: one just multiplies the autocorrelation function by the Fourier transform of ^ and adds the result to the original autocorrelation function. 5 APPLICATIONS AND DISCUSSION When discussing the application of the results we have ob- tained, it is important to interpret the results of observations correctly. Within the model we have used, all real observa- tional programs have two fundamental limitations. Firstly, it is the finite temporal resolution of observa- tions T, making the information on variability at high fre- quencies inaccessible to observations. For the power spectra, this can be modelled by multiplying them by some window function of the frequency which decays rapidly for u) > t~^ and does not affect it for low frequencies. A good choice here could be a Gaussian normalized such that its value at zero is unity: Wt{lli) = exp (53) Another important limitation is the finite extent T of any real observational program. This brings about an inability to estimate variations on very long scales and results in a power spectrum where adjacent frequencies are no longer independent; as a result, fine spectral features are smeared out. This can be modelled by convolving the power spec- trum with a similar window function of width ~ T~^; since the total power should stay constant, the function must be normalized such that its integral is unity. An example of such a function is an integral-normalized Gaussian WTii^,l-^ ) T 2-K exp T^iiu'-ijy (54) Let us now turn to some specific examples of the ob- jects and study how their variability is affected by gravita- 8 Tuntsov & Lewis tional microlensing. We will first consider periodic or nearly periodic sources, sketching the kind of effect gravitational microlensing brings about for objects observed through the microlensing population of the Galaxy on distances of order 1 kpc; it will be convenient to use the frequency domain in this study. Then, aperiodic sources will be considered with a spe- cific application to the extra scatter in the flux of stan- dard candles; here the language of autocorrelation functions will be more appropriate. An important point regarding the standard candles is that they are most interesting on the cosmological scales while the analysis present in this work is essentially Newtonian in its use of space and time. In ad- dition, the parameters of the lensing population on these scales can only be guessed; therefore the results obtained in cosmological applications should be viewed cautiously but do give a general feeling of what the effect of microlensing on these scales should be. In both cases optical (and shorter) wavelengths will be assumed, because at radio wavele ngths interstellar scintilla- tions are much more effective (e.g.. lScalo fc Elmegreerl2004i FLazio et al. 20041. 5.1 Pulsars The most prominent features in the power spectra are dis- played by objects which are periodic or close to being pe- riodic. One of the most interesting for us will be pulsars. The way their variability is affected is an interesting ques- tion on its own; in addition, they will serve as a tool to help us interpret the action of spectral broadening on this simple example. The power spectra of periodic objects such as pulsars display an array of spectral lines at multiples of the spin fre- quency ojp; the relative heights of these lines depend on the pulse shape. Broadening of every line is a typical feature of the multiplicative correction with a slowly- varying function. Indeed, let us consider a periodic source of flux J^p(i) = ^yl„e'""''*, A„ = yl! n with a power spectrum Ppiuj) = 2^ \^n\ S{u) ~ ncjp), (55) (56) modulated by some function H-e(t), with e varying so slowly that its Fourier transform e(Aaj) is essentially contained in the region Auj <^ uip. Then, for the power spectrum of Fp{t) — Fp{t){l + e(i)) one has^ Ppiio) = Pp{Lo) + Y^ \A„f P,{lo - nujp n + 2_^ A„A^e{uj — nu}p)e*(i (57) LU — mUJr, n^m where Pe is the power spectrum of e: Pe(Aaj) = |e(ALi;)|'^. The last line in l|57|l can be neglected due to assumption ^ We drop the infinitesimal correction to the discrete part of the spectrum. above, and comparing what is left to the results of our cal- culation in the previous section, we can interpret ^(Ao;) as the power spectrum of the magnification factor. Let us estimate the quantities involved for a typical pulsar at a distance of 1 kpc observed through the pop- ulation of white dwarfs of mass 0.5 Mq with the space density of 0.01 pc~'^ (as recent studies along different lines suggest - HolbcrE, Oswah & Sion '2002"; Hafizi ot al. '200^; .Qppenheimer et al. 2001 : Han fc Gould 2003') : assuming they have halo kinematics, we may choose o ~ 250km/s (e.g., iBinnev fc Trema ind Il98 q). For such w hite dwarfs, on/m » 3 X 10^ ijShapiro fc Texikolskvlll98,'Jl . Given that pulsars generally have high spatial velocities, a choice of v/a — 2 seems to be reasonable. Then, for InA = 30, the amplitude parameter becomes g « 17 and the frequency scale Aluq = 0.05 s~^. Given a relatively high value of the velocity parameter, tp starts off as an increasing function of Ao; as the latter changes from zero to Aluc ~ 0.01 s~^, where it turns to a shoulder followed by the y~^ tail for Alo > 0.3 s~^; a small step is present at AtjRi0.5s-^ The total extra power added to a line by gravitational microlensing is, according to 14911 . y « 15% which is con- tributed mostly by the shoulder - the central part accounts for about a tenth of this power while the tail makes up a fraction of order few percent. The tail contributes slightly more to the integral 1521 1. and the effective line width, as given by ||MJ, is (Acu) ^ 0.06 s~\ These frequencies scale according to 1G14II with param- eter a - e.g., the frequency difference scale would be twice as large if the defiectors' velocity dispersion was only a half of the one assumed and twice as low if the defiector population was two times faster. Our predictions for the effect of gravitational microlens- ing can be interpreted in a simple way for pulsars with spin frequencies greater than about ls~^. Namely, we suggest that the pulse height is modulated by a slowly (and, pre- sumably, randomly) varying function with an rms of 0.08 stellar magnitudes on the time-scales T > Auj~^ ~ 100 s; the time resolution required r < 1 s. The actual pulsars show a degree of pulse-to-pulse vari- ability that is much stronger than the expected 0.08™ con- tribution from lensing. However, unlike the lensing signal, these variations seem to be uncorrelated even on interpulse scale and therefore can be suppressed by averaging a large number of pulses. This does not necessarily mean a long time-scale - the lensing signal correlation time is of order seconds. Therefore, ceteris paribus, the lensing signal is more likely to be detected in high-frequency pulsars. 5.2 Quasi-periodic and related oscillations Another phenomenon of considerable interest in relation to the subject of this paper is quasi-periodic oscillations (QPOs) generally defined as broad - up to of order the cen- tral frequency - lines seen in the power spec tra of various astrophysical sources (e.g.. Ivan der Klial2005^ : more coher- ent features, with higher centroid frequency to width ratios Q, are known by a variety of names, deflned phenomenolog- ically and often specific to a class of sources or even to the particular state of a given source. Within the Galaxy, these oscillations are associated Microlensing in phase space II: Correlations analysis 9 with accretion onto compact objects - black holes, neu- tron stars and white dwarfs - from low-mass normal com- ponents. The frequencies of these oscillations span a wide range from mHz for cataclysmic variable stars that con- tain accreting white dwarfs to Hz and even kHz for low- mass X-ray binaries where the accretion proceeds onto neu- tron stars or black holes. Normally, features with Q in the range of up to a hundred or so are called QPOs, while for oscillations with greater coherence values terms such as normal, horizontal or flare branch oscill ations (N/H/FBOl in the case of low-mass X-ray binaries (Ivan der Klial2005i iHasinger fc van der Klij[l989l) and dwarf nova oscillations (DNOs), including longer-period ones (IpDNOs), in the case of cataclysm i c variables are in use IIKluzniak et alj l2005i IWarnejEooi : IWarner fc WoudlJl2004) . For convenience, we use the short term QPO more generally to mean all of the above when necessary. It is believed that careful modelling of these spectral features can hold potential clues to the physical processes in the inner part of the accretion discs and that QPOs may become important probes of the strong gravitational and/or magnetic flelds of the centra l compact objects (Psaltis 2001; iHeUieil Eooi : iMauchd l20o3) . Therefore, one has to inter- pret these features properly, including accounting for pos- sible contamination of the signal en route to the observer; obviously, gravitational microlensing, in the form described above, is one such possible 'contaminant' able to signifl- cantly modify the intrinsic power spectrum - although it does not change the the centroid frequencies of QPOs, the total power in the features and their width can be consider- ably affected by lensing. When making an estimate of the effect of lensing on the power spectra of these sources, the only change that seems necessary for the parameters of the lensing problem given in the previous subsection concerns the velocity - a choice of v/a = 1 is probably more sensible in this case. As a result, the central part is no longer suppressed but otherwise tp looks nearly the same. The total extra power added by lensing is slightly greater at y ^ 0.25, and the effective line width is lower: (Alo) ~ 0.035 s~^. We can therefore conclude that gravitational microlens- ing has the potential to influence the total power observed in a given feature of the power spectrum and should be taken into account in their interpretation. With regards to coher- ence estimates, lensing makes a negligible contribution to the line width of kHz QPOs but cannot be ignored in the interpretation of spectral lines at lower frequencies and may be the dominant effect in broadening lines at mHz frequen- cies. These results can be summarized in a graph that loosely classifies oscillations of different kinds in frequency- coherence coordinates (Figure [IJ. Apparently, whenever the lensing effective line width is of order 2-kv/Q or greater and y ^ 0{1), it should be taken into account. Therefore, lens- ing becomes important for Q > 27ri^/{Atj) - the regions that lie above the straight lines in the plot. 5.3 Aperiodic sources and standard candles Sources that do not show any prominent features in their power spectra, are more conveniently discussed in terms of the autocorrelation function. As mentioned in the previous Figure 4. A very rough classification of different spectral lines in terms of the central frequency u and coherence parameter Q = u/Au (Ai/ is the spectral line width). The left part of the plot is occupied by different features seen in the power spectra of cataclysmic variables while on the right one can sec more rapid os- cillations typical for accretion onto relativistic components. Grav- itational lensing becomes important in the interpretation of co- herence when broadening line width {Aui) becomes comparable to u/Q - these regions lie above the straight lines seen in the figure which are plotted for g = 5 X 10^, q = 17 and q = 0.5 correspond- ing to the distances of 10 kpc, 1 kpc and 100 pc, respectively; the other parameters are fixed as in the text. The total extra power corresponding to these three values of q are 0.32, 0.25 and 0.01 - that is, broadening contributes a little in the latter case although the line is formally quite broad. section, tf}{uj,uj') only depends on the frequency difference Alu = ij' — u! in the approximation we use. As a result, the correction Af{t) to the autocorrelation function f{t) can be found by multiplying f{t) by the Fourier transform of ^: A/(t) = /(t)V'W- (58) However, we should take into account the fact men- tioned in the beginning of this section - namely, limited temporal resolution of observations and thus the inability to observe high frequency signal. To model this, we multi- ply ■tl){A(jj) by the function Wt{i^) defined by (I53I I and then Fourier-transform the result to obtain an 'observationally- corrected' version of ip (we include \/27r factor in the defini- tion of ip): ■tPr{t) = / dije-'"^'Wr{co)i){uj). (59) Both tp and Wt are even functions oi uj, so the exponent in the last expression can be replaced with the cosine function. Denoting to = Aug^ we can rewrite the above formula in a form, similar to 14411 : Mt) - i^^-PJ-f/^-^l^, / * D In A to' Ds+D' a' with *r (?7,s,w,g) duj ■^{y,s,v,q)e yri-y^T^AL^i^/2 (60) (61) This integral also has to be done numerically. We have found the following rearrangement of the integral to pro- 10 Tuntsov & Lewis " ^^ S=0 v/a=1 10"' ^v ^ ^v \ A z=0.1 10"^ 10"^ S. \ ""■•-.. ^- 10- \ 10"^ \ \ — zxAco^= ZXACOq= =0 =1 \^ 10 ^ ,_,_ zxAco^= zxAco^= =10 =100 \ Figure 5. The autocorrelation correction ^^ as a function of rj = t/to for q = 10, s = 0, tj/cr = 1 and a few values of the time resolution parameter r/qto — zero, one, ten and a hundred. It has a flat central part followed by ri~^ tail. Figure 6. The height of the central plateau ^T\ri=o as a function of the time resolution parameter r/to for q equal to one, 0.1, 0.01 and 0.001; the velocity parameter ti/cr = 1 and s = 0. The quan- tity is weakly sensitive to the time resolution for t < (Aa;)~^ ^ toq/lnq and decreases oct^^ after that. duce fast and stable convergence for various values of the parameters (again, we assume A is so large that the Bessel function it appears in can be neglected): *^(r?,s,ii,g) (62) dx X wi^s, x) Jq (x/q) j{rjx/q, rAujQx/q, v) with an auxiliary function j{h,t,v) dy cos hy exp 1 + i 2 — ?; — y loivy). (63) At u = 0, q is simply a Gaussian of width \/T+t? while at non-zero v, this function displays a few oscillations around zero value before being exponentially suppressed on the same scale ~ O (\/l + t^j . Based on the similarity between j{h) and a Gaussian, one can expect that ^t has an effective width of order (Aa;)~^, which is approximately ~ toq/lnq for g > 1. Simi- larly, we expect that r should assume the role of this char- acteristic time when it exceeds the latter; cutting off some power at higher frequencies also leads to certain suppression in the value of ^t at low t. Figure confirms these expectations. It shows the be- haviour of ^I/t as a function of r; = t/to for q — 10 and a few values of the parameter r/io - 0, 10, 10^ and 10"^ that correspond to t{Alo) of zero, approximately one, ten and a hundred. One can see a flat central plateau extending to approximately the greater of (Aa;)~^ and r, followed by a ri~^ tail which is nearly independent of r. The height of the plateau as a function of r/to for a number of q values is shown in Figure 1^1 In Figure |51 a small value of the velocity parameter v/a = 1 is assumed. For higher velocities the Fourier integral for $ is dominated by non-zero frequency modes, and ^t drops below zero outside the flat central part, then tends to zero value in the same, rj Figure |7| illustrates this. ', manner from the negative side. \ tAM|j/q=1 \ tAm /q=3 ,,, TA(O|j/q=10 1 1 1 1 1 I \ s=0 ^ 1 v/<j=5 . I 1 t 1 \ \ q=10 ^ l/C-- ,' Figure 7. Behaviour of the autocorrelation correction <&,- as a function of r] for a high value of the velocity parameter v/a = 5 and a number of different values of the time resolution parameter r/qto — one, three and ten; the amplitude parameter g = 10 and s = 0. To incorporate both the different scales and the negatively- valued region, the linear scale is used but the functions are divided by their values at r; = 0. Although not easy, it is interesting to compare these findings to the results of previous works on calculation of the lensing au t ocorre lation function. Both Seitz et al. (1994) and iNeindor J (|200.'l ) perform these calculations for a two- dimensional distribution of lenses, and the resulting curve depends on a number of parameters, including the effective convergence and shear and scaled radius of the source, while this work assumes a (nearly) three-dimensional distribution of lenses with no external shear, and calculations are per- formed for a point-like source. However, inspection of Figures 9 and 11 of Seitz et al. (1994) shows that their autocorrelation function r{Y) does drop off as approximately an inverse square of Y in the shearless and parallel-to-shear cases for Y > 3 while Microlensing in phase space 11: Correlations analysis 11 perpendicular-to-shear curve is clearly not a power law; it should be noted, however, that this impression is based on measurements of the actually published plots. On the con- trar y, the growth of the autocorrelation function in Figure 1 of lNeindora i2003r seems to be shallower than the inverse squares; however, it is the result of a best-fit procedure to the actual autocorrelation function observed in one of the quasars. Let us now describe how an autocorrelation function of a distant source is affected by the transformation we have presented. For any source, the intrinsic flux autocorrelation function fit) should tend to the square of the flux average value for large i. Let us for simplicity assume this value to be unity: fit) —* 1 for t — ^ cxd. On top of that, one has all the information about second-order statistics of variability, normally presenting itself as a decreasing - not necessar- ily monotonically - function of t. In particular, the t — value of the autocorrelation function gives the flux disper- sion: cr^ = /(O) - 1. When / is subjected to the transformation described in this paper, the correction is proportional to /(i), not the variable part /(t) = f{t) — 1, and as a result, even the con- stant brightness sources are affected: fit) = fit) - 1 = m + m (1 + m) ■ (64) As can be seen from 16011 and the results of the calcula- tions in this section (e.g.. Figure |HJ, t/j is less than unity for any realistic In A. Therefore, for mildly variable sources with /(O) = cr'^ < 1, the above formula can be rewritten in the most simple form fit)^nt)+m (65) - that is, the microlensing correction simply adds up to the intrinsic autocorrelation function. The amount added is actually rather modest. For In A = 30, the factor in front of exp[— u^/2o-'^]^ in 1601 is approx- imately 0.41 and ^(0) saturates at about 1.3 for z <^ 1. As a result, ijj{t) < 0.5, and as with the intrinsic fit), the maximum is at t = 0. For sources with some memory of their past, one expects fit) to be smooth at f = and therefore the '(/'(t) plateau is hard to spot. Whether gravitational lensing becomes a dominant contributor to /'(i) therefore depends largely on the shape of the intrinsic function fit). As ip cc t~^ at large t, it necessarily becomes dominant if fit) decreases as a power law f~" with n > 2 or more rapidly - e.g., exponen- tially. The well-sampled observed autocorrelation functions are indeed often fit with exponentials but t his is mostly for the purpose of simp licity and convenience llHawkinj|l996i Ide Vries et alj|2004f) . A simple transformation of the auto- correlation function - the structure function - is more often used to characterize variability, but despite its simplicity, it is non-linear on /(f) and therefore is less suitable for our analysis. The power spectra of various quasars and AGNs have also been studied o bservationally in a number of works fGiveon et al." '1999'; 'P altani et all Il997l : iFiore et alj Il998l: [Edelson fc Kirpal 19951) . Traditionally, their variability is divided into a short and long time-scale; the former has stronger amplitudes and is thought to arise from the pro- cesses in the inner engines of quasars while the origin of the long-term variability is less clear and a variety of differ- ent mechanisms have been proposed (e.g.. lBregmanlll99d : IZackrisson et al.ll2003r : it is quite likely, that a number of mechanisms actually contribute and microlensing could be just one of those. We therefore are only concerned with the long-term variability. According to the theory presented in this study, any strong signature in the power spectrum of a quasar will be broadened by a universal function i/)(Aaj). In particular, if a quasar or AGN has only one dominating line at tj = 0, its entire power spectrum will be given by this function. According to estimates given in the previous section, for q ^ I most of the power resides in the Auj~^ shoulder of tp. The observed power spectra on long time-scales are well described by a featureless power law with and index ~ —2 outside the central region where the spectrum flattens (it is hard to estimate it here, though, due to the limited ex- tent of the light curves) . Edelson & Kiroal ( 1999|) give an estimate of —1.74 ± 0.12 for the slope of the X-ray fluctua- tions power spectrum of a nearby Seyfert galaxy NGC 3516; the characteristic frequency of the break in the spectru m is ujc « 27r/l month ^ 3 x 10"*^ s~\ Iciveon et alJ lll999ll re- port an average slope of —2.0 ± 0.6 (at frequencies range 10"^ s~'^ <uj < 10"'' s"^) for the sample of ~ 40 Palomar- Green quasars observed regularly in the optical over seven years; based on the simulated sensitivity of their method, they conclude that they cannot reliably infer these quanti- ties for individual quasars in the sample. They do, however, report two cases with an apparent break in the power-law slope at a frequency cuc ~ Sx 10~^ s^^. However, estimation of the power law index is still difficult with the present-day data accuracy and coverage. The total variability power in these spectra are of order a few to about twenty percent (rms). Therefore, although it would be extremely preliminary to claim that these data support our predictions, they certainly do not rule out the microlensing hypothesis at the present time. Some other aspects of this variability, such as chromaticity and (tentative) asymmetry of light curves, however, arc more difficult to accomodate within it (Zackrisson et al. 2003). Nevertheless, it would be interest- ing to consider whether the microlensing-induced broaden- ing, if ever detected, can place meaningful constraints on the global microlensing population. The model presented in the previous section is probably too crude to describe the actual distribution of microlensing bodies in a cosmological situation with strong clustering of lenses, and more work needs to be done to infer quantitative characteristics of this population from the available data. Even if microlensing never shows up as a dominant com- ponent of the autocorrelation function or the power spec- trum of a source, it still contributes significantly to the dis- persion of the flux. Indeed, *I'(0) saturates very quickly with g > 1 (see Figures 13 and |SJ and for the lensing population considered in this section, this happens relatively quickly, at distances of order a few kpc. As a result, for moderate values of v/a, one gets a top-up to the magnitude of order 0.1™ — 0.2™, as estimated in subsection 5.1. In the language of frequencies, it is clear that this vari- ability is spread over a range of timescales from roughly to to qto and consequently can be detected only if the time resolution r of observations is better than to; a fraction of order r/iqto) of this variability will be 'smoothed out' if 12 Tuntsov & Lewis r > to. The larger time-scale T is responsible for the sample completeness - for shorter light curves large variations are likely to be left outside its coverage; this is not an issue if the dispersion of an ensemble of similar sources is consid- ered, which is the case when we calculate the extra scatter in the brightness of standard candles due to microlensing is. Perhaps, the most interesting of such questions deals currently with the magnitude scatter of distant type la su- pernovae used in cosmology to measure the luminosity dis- tance as a function of redshift for constraining the evo- lution of the Universe JRiess et alJIlOQa : IPerlmutter et alJ Il999r . The full lensing aspect of th e problem is beyond the scope of this paper (see, e.g.. iHolz fc L indcr 2005; lAmanuUah. Mortsell fc GoobaiJIJOol iHui fc Green..2005., for recent studies of the role lensing plays here) as the geometry of the Universe assumed here is flat and static and cannot describe the high-redshift Universe. However, the treatment of this paper still applies to the Universe at moderate redshifts up to at least 0.01 - 0.1, and this is precisely the redshift range used to determine the zero-point of the rcdshift-magnitude relation and its uncer- tainty ( Wood-Vascy ct al. 2004) ; an estimate of the effect in cosmological settings should also give some idea of what kind of time-scales are involved. The cosmological microlensing population is still quite a mystery. For the economy of the hypotheses we will still assume that the deflectors are unseen white dwarfs of 0.5Mq which are uniformly distributed with density v. It can be written in terms of the fraction of the critical density in deflectors fid and Hubble constant Ho as: _ SHo^d which leads to 2Dy/Q.d\n\\/mD Dt Po (66) (67) where the Hubble distance Dh = c/Ho has been introduced; one can also recognize yrnD as the radius of the Einstein- Ch wolson ring for a lens of mass m seen from a distance D llSchneider. Ehlers fc Faicolll992l) . D/Dh is very close to the redshift z for small values of the latter in a flat Uni- verse. For a c urrently preferred value Dh JLahav fc Lid dlc 2004), one has ^"'"^ %o/10'^km/s M Mq 1/2 3/2 (68) An estimate fid In A ~ 1 can be regarded conserva- tive since we know that MACHOs contribute around ten to twenty per cent of t he dark matter in ou r neighbourhood jAlcock et alJ |200C| iLasserre et alJ l2000l : ISadouled Il999t lEvans fc Belokurov^l2004^ ■ and the dark m atter share in the comp osition of the Universe is Q,dm ~ 0.27 JLahav fc Liddld 12004^ . As a result, g ^ 1 at redshifts z ~ 10~^ and above - which includes all intergalactic scales. It also seems sensible to choose the same velocity dispersion a = 250km/s (e.g., iTormen et al]ll993l: iBahcaU fc Siandll996ll and the velocity parameter v/a = 1. For these parameters, V'(O) _„ ~ y ~ 0.35. The time- scales are to = ^'-^n = 20 s and gio ~ 8 X 10*^ s Vf^dlnA po/lO'J cm M Mq 1/2 3/2 (69) At z = 0.01 the latter time-scale is approximately 10 days. The brightness of a supernova is normally derived from a fit to the obs erved light curve with the extent of order ten to fifty days JRiess et alJl2004r . which, therefore, should be considered the time resolution. Given that Vt(0) oc r~^ at these scales, we can expect the extra scatter in the magni- tude to be MO) 2 -0.4 In 10 0.04" 0.2" (70) at these redshifts. Such an extra scatter is comparable to the intrinsic scatter in the brightness of type la supernovae, which is currently estim ated at am ~ 0.15"* iJTonrv et alJ 1200,4 IAs et al.ll200Rll . and therefore should be taken into account. At greater z ~ 0.05, which is roughly the target redshif t of the SuperNova Factory project ijWood-Vasev et al.l2004^ ■ the larger timescale becomes toq ~ 10^ days and the whole extent of the high-quality light curve typically used in esti- mating its brightness t < 50 days is not sufficient to signifl- cally alter the scatter due to gravitational lensing. Therefore, this project should be able to put some constraints on the spatial density of lenses v. In these settings, we cannot ex- pect these constraints to be particularly useful because most of crjji"'^ comes from the first ~ 10 — 100 kpc of the optical path. If a method could be devised to actually estimate the autocorrelation function from these data, it has the potential to serve such a probe. However, as mentioned above, other more stable objects seen to high redshifts (e.g., quasars) could be better targets for such studies. The numerical pre- dictions for the amplitude of the autocorrelation function correction and the time-scales involved are, of course, the same as for the supernovae as long as both can be consid- ered point-like. 5.4 Conclusions In the present paper we have calculated the autocorrelation function of the flux of a point-like source through a pop- ulation of uniformly distributed point-like gravitational mi- crolenses using the lowest nontrivial order approximation for the deflection potential defined by 11311 . In the last section we have also identified several potential applications in illus- tration of the method. Although, given the simplicity of the approximation, the estimates we obtained should be taken rather cautiously, they do show that the effect is important and cannot be ignored. It is sensible to conclude by recapping briefly the major assumptions used in this study, which should give us the idea of its applicability and show how and in which directions it can be improved. The 'physical' approximations employed here are the use of the geometric optics and thin lens ap- proximation as well as neglecting the potential time delays when writing down the propagation equation 0. The first of these was shown at the very beginning to be an excel- lent one in most situations in the gravitational lensing (e.g., Schneider et al. 1992), the other two were extensively dis- cussed in Paper I and Section 2 of the present paper. The major mathematical approximations are the trun- Microlensing in phase space 11: Correlations analysis 13 cation of Taylor series in 11311 on the first non-trivial, corre- lation, term and the perturbation technique used in the solu- tion of the propagation equation. The justification for series truncation comes from the fact that coefficients in front of the higher-order powers are all proportional to the appro- priate powers of the characteristic deflection angles, which only become of order unity very close to the lens; this never occurs for the deflectors considered in this paper; however, if one wants to consider black holes as deflectors, it will be nec- essary to look at these matters more carefully. The use of the perturbation technique can only be justifled retrospectively by comparing corrections to the uncorrelated propagation of the quantities involved. Here, T^ indeed turned out to be not more than about 0.5 of the zeroth-order quantity. If, for a difi^erent problem, this is not the case, one would have to sum the entire series 13211 or seek other solution methods. We have also only considered a point-like source and constant density of the deflectors. It does not pose any con- ceptual difficulties to perform the same calculations for a different source profile and deflectors density. However, such calculations are likely to require considerable numerical ef- fort. Apart from these approximations, we have also made certain assumptions about different parameters of the lens- ing population and most of our numerical predictions would change if very different parameters were considered. How- ever, there are some generic features which do not suffer from these uncertainties, the most robust of them are the shape of the power spectrum broadening function ■!/'(Atj) and the autocorrelation function correction V'(i)- In addition, due to the saturation of the variation amplitude at relatively low q ~ 10, it is not very sensitive even to considerable changes in v, po or m. To sum up, we have shown that microlensing needs to be taken into account in the interpretation of the power spectra and variability amplitudes of various sources. These predic- tions are well substantiated in the case of the Galaxy, and therefore timing properties of pulsars and X-ray binaries has the potential to reveal properties of the microlensing pop- ulation of the Galaxy. For cosmological distances, our pre- dictions are weaker because the knowledge of the microlens- ing population is less certain and the model of uniformly distributed lenses is unlikely to be applicable in this case. When this difficulty is overcome, one can attempt to study the cosmological microlensing population via this method. Acknowledgements. AVT is supported by IPRS and IPA from the University of Sydney. We thank Rodrigo Gil- Merino for useful comments regarding the manuscript and the referee for suggestions which lead to a considerable clar- ification of the presentation. AVT also acknowledges owing drinks to J. Berian James for the effort he put into the work on the English of this paper. REFERENCES Abramowitz M., Stegun I.A., 1972, Handbook of mathe- matical functions, 10th edition, 503-536 Alcock C. et al., 2000, ApJ, 542, 281 AmanuUah R., Mortsell E., Goobar A., 2003, A&A, 397, 819 Astier P. et al., 2005, A&A accepted l |astro-ph/0510447| l BahcaU N.A., Siang P.O., 1996, ApJ, 462, L49 Binney J., Tremain S., 1988, Galactic Dynamics, Princeton University Press Bregman J.N., 1990, Astron. and Astrop.Rev., 2, 125 Chandrasekhar S., 1943, Rev. Mod. Phys., 15, 1 Deguchi S., Watson W.D., 1987, Phys. Rev. Letters, 59(24), 2814 de Vries W.H., Becker R.H., White R.L., Loomis C., 2005, AJ, 129, 615 Edelson R., Kirpal N., 1999, ApJ, 514, 682 Evans N.W., Belokurov V.A., 2004, preprint (astro-ph/04 11222) Fiore P., Laor A., Elvis M., Nicastro P., Giallonco E., 1998, ApJ, 503, 607 Giveon U., Maoz D., Kaspi Sh., Netzer H., Smith P.S., 1999, MNRAS, 306, 637 Gradshteyn I.S., Ryzhik I.M., 1994, Table of integrals, sums and products, 5th edition Hafizi M., De Paolis F., Ingrosso G., Nucita A.A., 2004, Int.J.Mod.Phys.D 13, 1831 Han C., Gould A., 2003, ApJ, 592 172 Hasinger G., van der Klis M., 1989, A&A, 225, 79 Hawkins M.R.S., 1996, MNRAS, 278, 787 Hellier C., 2004, in Tovmassian G., Sion E.M., eds. Com- pact Binaries in the Galaxy and Beyond, Rev. Mex. de Astr., Vol. 20. Holberg J.B., Oswalt T.D., Sion E.M., 2002, ApJ, 571, 512 Holz D.E., Linder E.V., 2005, ApJ, 631, 678 Hui L., Green P.B., 2005, preprint (astro-ph/0512159) Katz N., Balbus S., Paczyhski B., 1986, ApJ, 306, 2 Kayser R., Refsdal S., StabeU R., 1986, A&A, 166, 36 Kluzniak W., Lasota J. -P., Abramowicz M. A., Warner B., 2005, A&A, 440, L25 Lahav O., Liddle A.R., Eidelman S., 2004, Phys..Lett. B 592, 1 Lasserre T. et al., 2000, A&A, 355, L39 Lazio T.J.W., Cordes J.M. de Bruyn A.G., Macquart J.-P., 2004, New Astron. Rec, 48, 1439 Lewis G.F., Miralda-Escude J., Richardson D.C., Wambs- ganss J., 1993, MNRAS, 261, 647 Mauche C.W., 2002, ApJ, 580, 423 Neindorf B., 2003, A & A, 404, 83 Oppenheimer B.R., Hambly N.C., Digby A.P., Hodgkin S.T., Saumon D., 2001, Science, 292, 698 Paczyiiski B., 1986, ApJ,304, 1 Paltani S., Courvoisier T.J.-L., Blecha A., Bratschi P., 1997, A&A, 327, 539 Perlmutter S. et al., 1999, ApJ, 517, 565 Psaltis D., 2001, Adv. in Sp. Res., 28, 481 Riess A.G. et al., 1998, AJ, 116, 1009 Riess A.G. et al., 2004, ApJ, 607, 665 Sadoulet B., 1999, Rev. Mod. Phys., 71, 197 Scalo J., Elmegreen B.J., 2004, Ann. Rev. A&A, 42, 275 Schneider P., Ehlers J., Falco E.E., 1992, Gravitational lenses. Springer Schneider P., Weiss A., 1987, A&A, 164, 237 Seitz C., Schneider P., 1994, A&A, 288, 1 Seitz C., Wambsganss J., Schneider P., 1994, A&A, 288, 19 Shapiro S.L., Teukolsky S.A., 1983, Black Holes, White Dwarfs and Neutron Stars, Wiley Tonry J.T. et al. 2003, ApJ, 594, 1 14 Tuntsov & Lewis Tormen G., Moscardini L., Lucchin F., Matarrese S., 1993, ApJ, 411, 16 Tuntsov A.V., Lewis G.F., 2006, MNRAS, this volume (Pa- per I) van der Klis M., 2006, in Lewin W., van der Klis M., eds. Compact Stellar X-ray Sources, Cambridge Univer- sity Press van der Klis M., 2005, Astron. Nachr. 326, 798 Wambsganss J., Witt H.J., Schneider P., 1992 A&A, 258, 591 Warner B., 2004, PASP, 116, 115 Warner B., Woudt P.A., 2005, in Hameury J.-M., Lasota J.- P. , eds. The Astrophysics of Cataclysmic Variables and Related Objects, ASP Conf. Series, Vol. 330. Witt H.J., 1993, ApJ, 402, 530 Wood-Vasey W.M. et al., 2004, New Astr.Rev., 48, 637 Zackrisson E., Bergvall N., Marquart T., Helbig P., 2003, A&A, 408, 17 APPENDIX A: UNPERTURBED BASIS FUNCTIONS AND INITIAL CONDITIONS A complete eigensystem of solutions to the propagation problem for one-ray statistics was found in Paper 1: exp [— i (k • r + ujt)] with M^Lm = c;™(t) (27r)3/2 Ak^nm = ^t;— k^ +i-yiK(cJ) (1 + \n\ +2m) . 2lu c where {x '■ ''' ~ t^x) , exp[— i(nx + ^k ■ r)] 27r and ^nm \'^ ) — /2m!\/iK(tj) (m+\n\y. 1 ( {/iK(a \n\ (Al) (A2) (A3) (A4) X exp y/iK{uj)l -\n\ (\/M^t-' The 'oscillator constant' n in the above expression is defined as Kic) ^ :^\\c\\{f3') (A5) and, according to our convention, when a square or quartic root of in is calculated, one takes the value in the complex plane closest to the positive ray of the real axis. The conju- gate basis Af^J^^ can be obtained by complex-conjugating Cnm(T) function only. See Paper I for details. All possible combinations of the form Ma{ri,r2,ti,t2,Ti,T2) = M^^^ (n , ti , Ti) M^^' (ra, ^2 , Ta) and associated sums \a = Ajj^ +^a2 with a composite index a = (ki,2, ^ =F '^/2, m.i^2,ni^2) form a dense set of solutions to the unperturbed problem 1251 1 with zero correlation term. In this basis, the initial conditions for a synchronous isotropic source of the form I12II can be written down as: {mIm^) 16nc^ g(ki)ff(k2)P(a;) cu'^D^, X exp - (ki - k?) D. (A6) Q-f 1 C(l ijjDs\J n({jS) (A7) APPENDIX B: FLUX FUNCTIONALS In order to obtain the expectation value of the flux second moment one needs to project the intensity moment vector 71/2 onto T = 0: F2 = M2|^^ ^^0 (Bl) As explained in the main text, this projection leaves intact the spatial and temporal sectors and therefore for ki,2 = K =p k/2, ijj\^2 = f2 =F t^/Z, Fourier modes of the flux are given in Ma basis as a sum F2(K,k,n,c^) = ^c4D)e^"^CaU,,,.o (B2) with respect to angular modes only. In particular, since C^ranij = 0) = S^d \JiK.{uj)/-K, one has for fi = 0: F2{K,k,Lu;D) = 71" Z — / mi. 2 CKk^mi.,(I?)e^'^''-'"i2-0(B3) Note, that ni,2 7^ can be ingnored. These moments can be analogously split into the zeroth- order contribution F2(K,k,(^;D) = } Ck TV ^ — ' CKko; (B4) and a term due to the correlated deflections (see II36II ') Fi{K,k,u;-D)= (B5) VM^) dD'Y^^e'^^'^-'^'^Ca'cl-e'-"'', where b = (K, k, w, 7711,2). The summation above is, in general, rather difficult to perform. However, significant simplification occurs for an isotropic and synchronous source of form (IA6II . First of all, initial conditions for isotropic sources are propor- tional to 5^,, hence one does not need to perform summa- tion in n' (these components do develop for D > but do not contribute to the fiux, according to JB3l l). Second, the set cjj in this case is completely defined by it spatial pro- file g(K' -'k'/2)g{K' + k'/2) and power spectrum P{uj'); hence, due to the linearity of the problem there should be coefficients Fad'{D) such that: F(K, \i,u;D)=J2 P^d' (^)^(^')fl(k'i)5(k': (B6) with d — (K,k,(jj); here we use the summation instead of the integration sign to avoid unnecessarily long notation. Noting that for sources described by IA6II e^'""!"^^ =e^'*^,7,(D)™i (773(1?))'"^ (B7) (B8) cLi^, = p('^)ff(ki)<7(k2)c«rHr)'^ with rid{D) = exp a-1 -7^(1 + i)D ^ c°.= a + 1 Microlensing in phase space 11: Correlations analysis 15 (B9) (BIO) and Ad i^(K.k) VM^: (Bll) (B12) it is easy to see that in the zeroth and first order the coeffi- cients Fdd' are given as Fdd' and Orf' TV ^ — ' 1 _ \/M^) -fdd' — Cjj/ J 7-i' ->>rf(D-D') A,/D' dD e ''^ 'e <* (B14) E/^ (^ l\vn-. / /-* /*\m^ 7n-\ / *\m.9 where r] = rid{D — D') and -q' — rid'{D') and the sum is over w-i,2,?Tii,2 running from to infinity. Summation in lB13t can be done immediateiy using X; z"' = 1/(1 - z) vaiid "iz e Z : \z\ < 1 and gives F'L'iD)^ (B15) d 16c ,2 exp -i2£ (K • k) {Ds +D)-^ V^^ a + 1 — (q — 1) exp ^^V^H(l + i)^ In order to evaiuate sum in lB14t we recali that Co is a function of r and t only and therefore the matrix elements Cw = -^||/:||(/3') (Aft, n^ Cora Mfc) (B16) a'^.t- t (^ \r l^ exp(iQi-Ti) X / d T2C+ij/2,0,m2(T"2)C+^'/2,0,m' (t2) 27r exp(iQ2 ■ T2) xr/ Co(k - k',a; - uj')t2 ( --||/:| |(/3")5(K - K'] where Qi T2c K =F k/2 K' =F k72 (B17) and C'o(Ak, Alj) is defined by 12111 . Functions ^u,o,m are given by (IA4II : Ca;,0,m(T) = CS-'^m f •\/iK(c<j)7 (B18) with the coefficient CZ = y 2^iK(Lj) exp — |y'iK(a7)r^ dependent of m. Therefore, one can change the order of summation with respect to m and integration with respect to r in 1B14I I and make use of the Laguerre polynomials generating function: 1-t exp xt 1-t E^'"( x)t" (B19) to obtain for the sum in the second line of IIB141 J(D') ^-'iu\\C\\{lf)^K(Lu)K{Lo')5{K-K!) (B20) exp [-7riV2 - -y'T^^] j2 j2 d Tid T2 (l-r,)(l-^»?')(l-r)(l-e'7'*) X (27r)-^c'«-^i'n^Co (k - k', a; - J) rae'^'^-^^^ where ^ , /IT-TTTi+i!/ 7 = (1 - i) vM^Y^ + vM^j -^r,' (B21) If one now denotes yii,2 the directions of Qi,2 <y5i,2 : Qi,2 = Qi,2e^i2, angular integration in l)B20^ yields a product of two first-order Bessel functions and radial integration in r can be performed using iden- tity fp°°drTVi(QT)exp(-7 rV2) = g/72 exp(-gV27) JGradshtevn fc Rvzhik^ll994^ . Thus we arrive at I{D') = AvyfK,{uj)K(Lu')5 (K - K') C (Ak, Alj) (B22) Q1Q2 exp [— Q?/27 — Q2/27*] cos(v3i — a) cos(q — 1^2) 72(1 - r,)(l - ^V) 7*2(1 _ ^.)(1 _ eV*) and , _ ,^<:Jr 16c^ exp[-ig£(K'.kOO.] -f'dd'(W) - \/ -7-77 -777^ n , ,.,12 (B23) fi:(a;') Lj'2i32 1 + a- I AD' I{D') exp [Xd{D - D') + Ad-D'] APPENDIX C: MODIFICATION OF THE TEMPORAL POWER SPECTRUM In this appendix we calculate quantities 13911 and the line profile function %j) in the simple case of a point-like source at the origin of the source plane (?(k) = 1 as seen by a single observer at the origin of its plane ri,2 = 0. Using IB15II one obtains, in this situation (care must be taken when calculating the limit of k — > because a tends to infinity in this limit): T\uj',D) 4P{lo') (CI) exp -fV^Ho Df {Ds + DY a + 1 — (q — 1) exp -fy^(l + i)Z) P{lo') {Ds + D) -^o{^). Therefore, in the leading order, the line profile function 1/) is given by the expression 1)41^ . The last line of lUH also shows that the zeroth-order correction to the free propagation of fiux is of the second order in k (and therefore v); as shown below, correction in the first order of perturbation due to correlated defiections is O {y/K\ and therefore represents the dominant term. To perform integration 14H for the first-order correc- tion, we first untilize the 5-function in K — K' from 1B23I I and then, taking into account that C depends on Ak = k— k' only, perform a change of integration variables (K, k, k') ^ 16 Tuntsov & Lewis (Qi, Q2, Ak) according to UJLO K k = k' 4cAlj LOLO (Q1-Q2) -^(Qi+Q.) + ^Ak (C2) (C3) The Jacobian determinant of this transformation £'(K,k,k') _ / uuj' ^ " D(Qi,Q2,Ak) ^ \^2cAcj is independent of integration variables and thus appears as a factor in front of the integrals, and since this is a linear transformation Qi,2,Ak independently run through entire VJ^ planes. Then, we may rewrite the first-order correction 14H as i:^{LO,J ,D) = k(cj) (C4) {Ds + DY ■ exp /a.- 4i^A/K(a;)K(a)')exp -D' {^^J k{uj') - ^a/k(cj)) (l-»?)(l-C'?')(l-r)(l-e'7'*) X d^AkC(Ak,Aa;) / dQi dQ QlQl ry2^*2 exp 27 27* X / d(y9l d(y92 .2c COs((^l — Of) COs(</32 — Ct) 47r2 1— (K ■ k) (D - 7:)') .2c X exp -1— {iS^- k.}[U - U ) -\— (K' ■ k') (Ds + D') here we deliberately postpone some obvious cancellations to keep the asymptotic behaviour clear and provide for lens population parameters dependence on D. Using IIC2II we rewrite the argument of the exponent in the last line in terms of Qi,2 and Ak, angular integration yields a product of two Bessel functions and the identity used in deriving lB22t brings the integral with respect to Qi,2 to the following form: ' p^Afc^ u(Afc) = /2Aa;^ \^ uju' ) (Ds+D)- with (|7|^ = Im^7 + Re^7) ■ exp \—p Re7 Afc 1 p={Ds + Df I 2Alj^ ^ ^ h 5Im7 + 5 Re 7 (C5) (C6) where c Function ii(Afc) is independent of its argument vector direction, and integration of (^(Ak, Ao;) with respect to the angular coordinate produces Infeld (modified Bessel) func- tion 27r/o(uALi;/cr^Afc). Therefore, one has for the integral with respect to Ak: d^AkC(Ak)M(Afc) dAfc exp 1 2Aa;^ UJUj' 2 \ cr2Afc2 2v2Tvm p e 2a a (Ds + Df + p{j + 7*)AA: (C8) -'"(^)h(T-') -'•(*?-' The exponential factor decreases rapidly at Afc tending to both zero and infinity thus providing necessary convergence. It is reasonable now to introduce some approximations. It was s hown in Paper I that in physically reasonable situa- tions ^yK{Lo) <C 1. We can use this fact and expand 77, ^, r;' etc. in its power series retaining only the few leading terms as necessary. First of all, we have for the denominator in I{D') f ^2K{Lj){^^{D-D')y +o{k) and \l-iri'\^=2K{Lu')(^{Ds + D'] and similarly 1 lu'''DU{lu') l + a|2 ^ 2c2 o(k) (C9) (CIO) (Cll) For 7 the first term in expansion is purely geometrical and imaginary while the first real term is proportional to «:: -2ic- {D-D')+Lo'{Ds ujuj'{D- D'){Ds + D') ^'ViO(^^) +ijJK[uj) Ds + D 3c D-^ {Ds + D){Ds + D') (C12) + 0(«'), where we have used u)k(u}) = u)'k,{u)') from IIA5II . Expanding Auj^ and using the above formulae one finds that at this level of approximation p is purely geometrical in the leading order: p = {D~D'f{Ds+D' (D, + i?)2 + 0(^2 {Ds+Dyr{i-ty+o{K){ci3) with dimensionless distance coordinate t = {Ds + D')/{Ds + D) introduced; as D' changes from Ds to D, the correspond- ing range for t is from s = Ds/{Ds -I- D) to 1. Now, assuming that parameters v, m, a, v are indepen- dent of D', we use HASH and {/3^) = 27r||/:||"^m^ In A, intro- duce the frequency difference scale Aujo 2a Po' dimensionless variables x = ^Afc, 2 Auj y Aluo and an amplitude parameter 2ljk{uj) 2 (Ds + Df'^ 3c Po (C14) (C15) (C16) (C17) m /87rlnAi/(Ds + -D) Po we write down the following approximation for the ra- tio llC4|l : ^{Auj,D)^ 12 exp [-uV2o-^] 2tt In A Aluo (C18) Microlensing in phase space 11: Correlations analysis 17 Figure CI. Auxiliary function w{s,£) defined by the inte- gral iC20i . Behaviour of «) as a function of J is displayed for five different values of s - 0, 0.2, 0.4, 0.6 and 0.8; w(s,i) oc ^"^ for large 1^. Xq I <it <ixt (1 — t) exp 2 2 -— -(1-i) {t -s ) X exp y 2x2 ^0 f-^) [Jo (x) - Jo {Ax)f + o{V^). - as expected, it is 0{y^). From this point, integration needs to be numerical, and is, perhaps, most easily done by first calculating the integral of the function on the second line of IClSt with respect to t and then performing calculation in x. This approach yields the expression 14411 with the scaled line profile function ^(y,s,^,q^ y '2x^ q I dxw {s,xq) exp defined through an auxiliary function 1 w{s,^)= dtt^{l~tfexp (C19) Jof-^) [Jo {x) - Jo {Ax)f \a X J ic'(l-i)'(i'-s')' (C20) The behaviour of w(s, ^) as a function of ^ is shown in Fig- ure lUll for a few values of s; w(s,5) oc ^~^ for ^ S> 1. Since A S> 1 we can safely neglect the second Bessel function in the square brackets when evaluating (IC19II .